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Abstract 

We establish universal behavior in temperature dependencies of some observables in 
(s + id)-wave BCS superconductivity in the presence of a weak s wave. There also could 
appear a second second-order phase transition. As temperature is lowered past the usual 
critical temperature T c , a less ordered superconducting phase is created in d wave, which 
changes to a more ordered phase in (s + id) wave at T c \ (< T c ). The presence of two 
phase transitions manifest in two jumps in specific heat at T c and T c \. The temperature 
dependencies of susceptibility, penetration depth, and thermal conductivity also confirm 
the new phase transition. 
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The Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity |], has been success- 
fully applied to different systems in pure angular momentum states such as s, p, and d waves. 
However, the unconventional high-T c superconductors [0] with a high critical temperature T c 
have a complicated lattice structure with extended and/or mixed symmetry for the order pa- 
rameter [ffl. Some of the high-T c materials have singlet c?-wave Cooper pairs and the order 
parameter has d x 2_ y 2 symmetry in two dimensions ||]. Recent measurements || of the pene- 
tration depth A(T) and superconducting specific heat at different temperatures T and related 
theoretical analysis |6|, [/]] also support this point of view. In some cases there is the signature 
of an extended s- or <i-wave symmetry. The possibility of a mixed (s — d)-wave symmetry 
was suggested sometime ago by Ruckenstein et al. and Kotliar |I0| . There are experimental 
evidences of mixed s- and rf-wave symmetry in compounds such as YBa 2 Cu 3 7 (YBCO) ||, 
and Bi 2 Sr 2 CaCu208+x @, where an [s + exp(i9)d] symmetry is applicable. Recently, this idea 
has been explored to explain the NMR data in the superconductor YBCO and the Josephson 



critical current observed in YBCO-SNS and YBCO-Pb junctions ||1 1|| . There have also been 



certain recent theoretical studies using mixed s- and <i-wave symmetries |L2| and it was noted 
that it is more likely to realize a stable mixed s + id state than a s + d state considering different 
couplings and lattice symmetries. 

It is quite natural that the Cooper electrons might interact in both s and d waves with 
different couplings. In the presence of simple central potentials the Cooper problem separates 
in its decoupled s- and <i-wave components. The same decoupling occurs in a linear Schrodinger 
equation. However, in the nonlinear BCS theory, the presence of both s- and ci-wave components 
in the interaction would lead to an order parameter of mixed symmetry and consequently a 
coupled set of BCS equations. The symmetry of the order parameter is to be specified in order 
to solve this coupled set of equations. 

The normal state of most high-T c materials has not been satisfactorily understood and there 
are controversies about the appropriate microscopic hamiltonian and pairing mechanism [||, [7] . 
Despite this, we study the (s-M(i)-symmetry case of the order parameter using the weak-coupling 
microscopic BCS theory based on the Fermi liquid model to extract some model-independent 
properties of such a description. For a weaker s-wave admixture, quite unexpectedly, we find 
another second-order phase transition at T = T cl < T c , where the superconducting phase 
changes from a pure <i-wave state for T > T cl to a mixed (s + z<i)-wave state for T < T cl . The 
specific heat exhibits two jumps at the transition points T = T cl and T = T c . The temperature 
dependencies of the superconducting specific heat, susceptibility, penetration depth and thermal 
conductivity change drastically at T = T& from power-law behavior (typical to d state with 
node(s) in the order parameter on the Fermi surface) for T > T cl to exponential behavior 
(typical to s state with no nodes) for T < T cl . The order parameter for the present (s + id) wave 
does not have a node on the Fermi surface for T < T c % and it behaves like a modified/extended 
s-wave one. The observables for the normal state are closer to the superconducting I = 2 state 
than to those for the superconducting I = state 0. Consequently, superconductivity in s 
wave is more pronounced than in d wave. Hence as temperature decreases the system passes 
from the normal state to a "less" superconducting ci-wave state at T = T c and then to a "more" 
superconducting extended s-wave state at T = T cl signaling a second phase transition. 

We consider a system of N superconducting electrons under the action of a purely attractive 
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two-electron potential in partial wave I (=0,2): 



V pi = ~ v i cos(W p ) cos(W q 



=0,2 



where 8 P is the angle of momentum vector p. 

Potential ([!]) for a arbitrary small V\ leads to Cooper pairing instability at zero tempera- 
ture in even (odd) angular momentum states for spin-singlet (triplet) state. The Cooper-pair 
problem for two electrons above the filled Fermi sea is given by [|7| 

vr l = £ co^ilO)^-^)- 1 (2) 

q(?>l) 

with the Cooper binding Ci = 2 — E\. Here e q = h 2 q 2 /2m with m the mass of an electron and 
the q-summation is evaluated according to 

. N r N r°° f 27T 
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Unless the units of the variables are explicitly mentioned, in this work all energy variables are 
expressed in units of E F , such that T = T/T F , = E^/Ep, E F = k B = 1, etc, where T F (E F ) 
is Fermi temperature (energy) and ks the Boltzmann constant. 

We consider a weak-coupling renormalized BCS model in two dimensions with (s + id) 
symmetry. At a finite T, one has the following BCS equation 

A P = -E^P^tanh^ (4) 

with E^ = {(e q — /i) 2 + | Aql 2 ] 1 / 2 . The order parameter A q has the following anisotropic form: 
A q = A + iA 2 cos(2#), where A/'s are dimensionless. The BCS gap is defined by A(T) = 
(Aq + A2/2) 1 / 2 , which is the root-mean-square average of A q on the Fermi surface. Using the 
above form of A q and potential (0), (f|) becomes the following coupled set of BCS equations 
for I = and 2 

1 = y cos 2(/m_L t anh% (5) 
Vi V 2E q 2T 

where the coupling is introduced through E^. 

Using @, set (||) of BCS equations can be explicitly written in terms of Cooper bindings 
as follows: 



dO cos 2 (Iff) 



2de r°° de E 

t 9 tanh q 



0. (6) 



2e q -E l Jo E n 2T_ 

The two terms in the BCS equation (|6]) have ultraviolet divergence. However, the difference 
between these two terms is finite. BCS model (|J) is independent of coupling VJ, is governed by 
Cooper binding C/, and has some advantages 0. Firstly, no energy cut-off is needed in this 
equation. This is why model (|6|) is called renormalized 0, [T3|]. Secondly, this model leads to 
an increased T c in the weak-coupling limit, appropriate for some high-T c materials 0. Here, 
we use the renormalized model for convenience. Otherwise, it has no effect on our conclusions 
and the same analysis can be performed in the standard BCS model with cut off. 
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The specific heat per particle is given by @ 



or) = ^ E / q (i - / q ) - (7) 



where / q = 1/(1 + exp(£' q /T)). The spin-susceptibility x is defined by 

x(T) = ^EA(i-A) (8) 



q 



where /ijv is the nuclear magneton. The penetration depth A is defined by 

2 



\-\T) = \- 2 (0) 



NT q 



EA(W q ) 



(9) 



The superconducting to normal thermal conductivity ratio K S {T) / K n {T) is defined by JF| 

ggg) S q fo ~ l)/ q (l ~ / q )gq 
K n (T) E q (e 9 - l) 2 / q (l " U) ' 



(10) 



We solved the coupled set of equations (H) numerically and calculated the gaps A and A 2 
at various temperatures for T < T c . The corresponding BCS gap A(T) was also calculated. 
For a very weak s-wave (d-wave) interaction the only possible solution corresponds to Ao = 
(A2 = 0). In order to have a coupling between s and d waves both the interaction potentials 
are to be reasonable. We have studied the solution only when a coupling between the two 
equations is allowed. In this domain we have kept the <i-wave coupling stronger than s-wave 
coupling, so that as temperature is lowered past T c a superconducting phase in d wave appears. 
In Fig. 1 we plot the temperature dependencies of different A's for the following two sets of 
s-d mixing corresponding to (1) C = 0.0006, C 2 = 0.001, (full line) and (2) C = 0.00085, 
C2 = 0.001 (dashed line), referred to as models sdl and sd2, respectively. For a superconductor 
with Tp = 5000 K, the largest of these Cooper bindings C 2 is 5 K. The smallness of this binding 
guarantees the weak-coupling limit, where the BCS model should provide a good description. 
In both cases the parameter A 2 is suppressed in the presence of a non-zero A . However, the 
BCS gap A(T) has the same form as in the case of pure s and d waves. In model sdl (sd2) 
A(0)/T c = 1.535 (1.644), T c = 0.0266 (0.0266), T cl = 0.01065 (0.0206). For a pure s (d) wave 
A(0)/T c = 1.764 (1.513) 0. At T = the order parameter has s- and <i-wave components 
and we find as T increases both components decrease and for T > T cl the s-wave component 
vanishes and one is left with a pure <i-wave component, which vanishes at T = T c . 

In order to substantiate the claim of the second phase transition at T = T cl , we study the 
temperature dependence of specific heat in some detail. The different specific heats are plotted 
in Fig. 2. With this two-step transition, the superconducting specific heat exhibits a very 
unexpected peculiar behavior. In both models the specific heat exhibits two jumps — one at 
T c and another at T cl . From (|7]) and Fig. 1 we see that the temperature derivative of |A q | 2 
has discontinuities at T c and T cl due to the vanishing of A 2 and A , respectively, responsible 
for the two jumps in specific heat. For T c > T > T cl , the specific heat exhibits typical ci-wave 
power-law behavior C s (T)/C n (T c ) = 2(T/T C ) 2 found in recent studies [0. For T < T c , we find 
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an exponential behavior. Two jumps in specific heat have been observed recently in certain 
superconducting compounds which suggest the existence of a coupled s + id phase Jf3| . 



Next we study the temperature dependencies of spin susceptibility, penetration depth, and 
thermal conductivity which we exhibit in Figs. 3 — 5 where we also plot the results for pure 
s and d waves from Ref. [|7| for comparison. In all cases ci-wave-type power-law behavior is 
obtained for T c > T > T cl . We obtain in d wave K S (T) « K n {T)(T /T c ) 1 - 2 and Xs(T)/ Xn (T c ) « 
(T/T C ) L3 0. For T < T c i, there is no node in the present order parameter on the Fermi surface 
and one has a typical extended s-state behavior. A passage from d to extended s state at T c \ 
represents an increase in order and hence an increase in superconductivity |7[] . As temperature 
decreases, the system passes from the normal state to a <i-wave state at T = T c and then to an 
extended s-wave state at T = T cl signaling a second phase transition. 

In conclusion, we have studied the (s+zd)-wave superconductivity employing a renormalized 
BCS model in two dimensions and confirmed a second-order phase transition at T = T c % in 
the presence of a weaker s wave. We have kept the s- and d-wave couplings in such a domain 
that a coupled (s + id)-w&ve solution is allowed. As temperature is lowered past the first 
critical temperature T c , a weaker (less ordered) superconducting phase is created in d wave, 
which changes to a stronger (more ordered) superconducting phase in (s + id) wave at T c \. 
The (s + id)-w&ve state is similar to an extended s-wave state with no node in the order 
parameter. The phase transition at T& is also marked by power-law (exponential) temperature 
dependencies of C(T),x(T), A\(T) and K(T) for T > T cl (< T cl ). A similar second-order 
phase transition may occur for some other types of mixtures of angular momentum states. 
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Desenvolvimento Cientifico e Tecnologico and Fundagao de Amparo a Pesquisa do Estado de 
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Figure Captions: 

1. The s- and d-wave parameters A , A 2 , and BCS gap A(T) at different temperatures for 
(s + ic/)-wave models sdl (full line) and sd2 (dashed line) described in the text with different 
mixtures of s and d waves. 

2. Specific heat ratio C(T)/C n (T c ) versus T/T c for models sdl (full line) and sd2 (dashed 
line). The dotted line represents the pure c?-wave result from Ref. |7| for comparison. 

3. Spin susceptibility ratio Xs(T)/x{T c ) versus T/T c for models sdl (full line) and sd2 
(dashed line). The dotted lines represent the pure s and d-w&ve results from Ref. |7| for 
comparison. 

4. Penetration depth ratio AA(T) = [A(T) — A(0)]/A(0) versus T/T c for models sdl (full 
line) and sd2 (dashed line). The dotted lines represent the pure s and <i-wave results from Ref. 
for comparison. 

5. Thermal conductivity ratio K S (T) / K n {T) versus T/T c for models sdl (full line) and sd2 
(dashed line). The short dashed lines represent the pure s and d-wave results from Ref. for 
comparison. 
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Figure 1 
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Figure 2 
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